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Abstract 



The stereoscopic reconstruction of air showers viewed by multiple imaging atmospheric 
Cherenkov telescopes (lACTs) allows a more precise reconstruction of shower energies and 
hence an improved determination of energy spectra. Reconstruction techniques and in 
particular new systematic checks are discussed on the basis of the large sample of 7-rays 
from Mkn 501 detected with the HEGRA lACT system. 
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1. Introduction 

One of the premises of stereoscopic imaging of air showers using multiple atmospheric Cherenkov 
telescopes is an improved reconstruction of shower energies and hence of source spectra. Stereo- 
scopic imaging allows a rather precise reconstruction of the core location; using the known 
relation between light yield, core distance, and shower energy, the energy is estimated by suit- 
ably averaging the measurements by the individual telescopes. The HEGRA system of imaging 
atmospheric Cherenkov telescopes (Daum et al., 1997) is operational since late 1996, and first 
results on energy spectra have been shown at various occasions; for example. Fig. ^ shows the 
current version of the Mkn 501 spectrum (Aharonian et al., 1997). Purpose of the present pa- 




Fig. 1: Typical energy spectra measured using the HEGRA telescope system. Errors are 
statistical only. 

per is to provide background information on the techniques used to reconstruct spectra, and to 

"'^Talk presented at the Workshop "Towards a Major Atmospheric Cherenkov Detector V", Kruger Park, 
South Africa, 1997 
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Fig. 2: Mean scaled width of events in different ranges of reconstructed energy. The full line 
shows the distribution of excess events from the direction of Mkn 501 (and a Gaussian fit), the 
dashed line indicates the off-source background. 

discuss the status of our understanding of the system and the systematic checks of the energy 
reconstruction, which are possible for the first time due to the redundant information provided 
by the telescope system. I will very briefly discuss the selection of 7-ray events and then con- 
centrate on two areas, the energy estimate and the determination of the effective detection 
area. I should emphasize that the primary goal in these analyses was to minimize systematic 
uncertainties, sometimes at the expense of some loss in statistical precision. 

2. 7-ray event selection 

For a 7-ray source such as Mkn 501, events are selected on the basis of their direction, of the 
event shape, and of the core location. The directional cuts and the estimate of their efficiency 
is relatively straight forward. Next, 7-rays are enriched using the shape information. The so 
called mean scaled width is calculated by determining for each telescope, as a function of image 
size, distance to the core and zenith angle the expected width of the image, normalizing the 
measured width to the expected width, and averaging over all telescopes. Fig. Elshows the mean 
scaled width for different ranges of the (reconstructed) shower energy. Excess events from the 
direction of Mkn 501 show up as a peak around a mean scaled width of 1; the cosmic-ray 
background gives rise to a broad distribution at higher width values. In the determination of 
spectra, typical cuts require a mean scaled width below 1.3, resulting in an efficiency well above 
90%. As an additional cut, the location of the shower core is restricted. Fig. IHl illustrates the 
distribution of shower cores as well as the locations of the four telescopes presently included 
in the lACT system. Event are required to lie within 200 m from the central telescope, CT3. 
Usually, we also remove events with a y coordinate greater than 150 m, reflecting the fact that 
the telescope CT2 is not yet included in the CT system. This selection guarantees that a least 
one telescope is within about 150 m from the shower core. 

3. Determination of core location 

A crucial input both for the event selection and the energy reconstruction is the reconstruction 
of the location of the shower core. The reconstruction is based on simple geometry (Daum et al.. 




Fig. 3: Distribution of reconstructed shower cores, together with the locations of the telescopes. 
Telescope CT2 (open circle) is not yet included in the lACT system. The dashed lines indicates 
the acceptance region. 

1997). With a 4-telescope system, it is not only possible to reconstruct shower cores, but one 
can also test the precision of the core reconstruction by dividing the full system into two groups 
of two telescopes each, reconstructing the core separately for each group and comparing the 
results. In Fig. El^a), the x-coordinate obtained from one group is plotted vs the x-coordinate 
from the other group. One sees a very clear correlation. Assuming that both measurements 
are uncorrelated, the error in the determination of the core location for such 4-telescope events 
can be determined (Fig. EJ^b)). The error in x (and also in y) is about 8 m for events with the 
core in the central region of the array, and increases to 20 m for distant events, where fewer 
telescopes trigger and where stereo angles are small. 

4. Energy determination 

As a first step in the reconstruction of shower energies, the response of the telescopes needs to be 
cahbrated. We use a two-step approach, first equalizing the relative response of the telescopes, 
and then adjusting the global scale. The determine the relative sensitivity of two telescopes i 
and j, shower cores are reconstructed using these two telescopes only. Then, the asymmetry in 
image sizes, (a^ — aj)/ {ai + aj) is plotted vs the asymmetry in the distances from the core to the 
telescopes, (r, — rj)/{ri + rj) (Fig. E)). If two telescopes have a similar sensitivity, the average 
asymmetry in image sizes is zero for rj vj (Fig. El^a)). In the telescopes are unequal, the 
distribution appears shifted (Fig. El^b)). We find the following ratios of sensitivities: CT4/CT3 
= 0.937 (0.937), CT5/CT3 = 0.857 (0.855), and CT6/CT3 = 0.791 (0.778). The first number 
is always determined using a cosmic-ray event sample, the second (in parenthesis) using the 
Mkn 501 7-ray sample. The two sets of numbers are consistent within their statistical errors, 
which range from 0.005 to 0.007. We believe, therefore, that the relative calibration is good to 
about 1% or better. The 20% differences between telescopes could be caused by differences in 
the quantum efficiencies of the different batches of PMTs, by variations in mirror reflectivity 
etc. (the telescopes came into operation sequentially, over a span of 1.5 years). 
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Fig. 4: (a) Scatter plot of x-coordinate of the core measured with one pair of telescopes, vs that 
measured with the other pair, for '-/-rays from Mkn 501 . (b ) Error in the determination of the 
x-coordinatc of the shower core for 4-telescope events as a function of core distance from the 
central telescope, derived assuming that the two 2-telescope measurements are independent. 
The resolution given refers to a Gaussian fit to the distributions; there are slight non-Gaussian 
tails. 
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Fig. 5: Asymmetry in measured amplitudes ai, aj for a pair of telescopes as a function of 
the asymmetry in the distances r^, rj from the telescopes to the shower core, (a) For CT3 and 
CT4, (b) for CT3 and CT6. 



The absolute calibration of the energy scale is much more challenging. Presently, we study 
three different methods: i) the determination of the photoelectron to ADC channels conversion 
factor using the width of the distribution of laser test pulses to estimate the mean number of 
photoelectrons per pulse, ii) the measurement of the response of the whole telescope and of 
the analysis chain simultaneously using a distant, pulsed, calibrated light source (Frass et al., 
1997), and iii) the calibration by comparing the measured and the expected cosmic-ray rate 
(A. Konopelko et al., 1996). For example, at a fixed reference wavelength of 430 nm the first 
method results in a sensitivity of (1.13 ± 0.25) • 10^ ADC channels/ (photon/cm^), the second 
in (1.20 ± 0.12) • 10^ ADC channels/ (photon/cm^). For the cosmic-ray comparison, the MC is 
still being refined, but first analyses indicate values consistent with the other two, with errors 
of about 15%. Ultimately, we believe that the energy scale can be determined within about 
10%. Note that all energy spectra shown so far use the energy scale derived from a previous 
generation of MC simulations, where the detector simulation had known imperfections; scales 
are likely to change by 10% to 20%. 

Once the sensitivity is known (or assumed), the shower energy is determined by converting 
the image size seen in each telescope into an energy value, using the measured distance to the 
core and the relation between energy, core distance, zenith angle and image size obtained from 
simulations. Two different approaches are used, one based on an analytical approximation, the 
other based on Monte-Carlo generated look-up tables. Fig.lHl^a) illustrates a typical dependence 
of average image size on core distance. To properly average over telescopes, the errors of the 
energy estimates need to be known. One contribution are shower fluctuations and photoelectron 
statistics (Fig. Efb)); these errors are large close to the core, and roughly constant at about 
20% for larger core distances. A second contribution arises from the error in the determination 
of the shower core; this contribution is proportional to the slope of the size-distance relation, 
and to the core error, and rises steeply for core distances of 100 m and above. 

A new and unique feature of the telescope systems is that these crucial inputs can be 
verified experimentally. To measure the light yield as a function of core distance, one selects 
events with cores at a fixed distance from telescope i, and with a fixed image size Cj. These 
conditions select showers with a (roughly) fixed energy. One then plots for one of the other 
telescopes j the image size as a function of the (varying) core distance rj. Of course, care has 
to be taken not to bias the results due to trigger conditions and selection criteria. Fig. [3 shows 
that the measured distributions agree reasonably with the simulations. 

The shower energy is derived by averaging over the measurements from the individual 
telescopes. Monte Carlo studies show that the resulting energy resolution is not sensitive the 
details of the averaging process, such as the exact choice of the weights, or if logE instead 
of E is averaged. Figs. |Hfa,b) shows energy resolutions slightly below 20% obtained with 
two different algorithms. Both show a non-Gaussian tail towards low reconstructed energies, 
caused by events where the core location is poorly determined. The tail, however, does not 
harm the reconstruction of steeply falling spectra. The energy resolution is almost independent 
of energy over a wide range (Fig. Efc)). Again for the first time the lACT system allows 
to check the energy reconstruction procedure by determining energies of 4-telescope events 
independently from two pairs of two telescopes each. Fig.^l^a) shows that the two measurements 
correlate quite well. Assuming that the two measurements are independent, one can derive the 
energy resolution of the system from the width of the distribution in the difference of the 
two measurements (Fig. El^b)). The resulting resolution figure of 9.3% is significantly better 
than the resolution expected on the basis of simulations, which seems to indicate that the two 




Fig. 6: (a) (Left) Typical relation between mean light yield and distance to the shower core, 
for TcV 7-ray showers, (b) (Right) Relative error of a single-telescope energy estimate as 
a function of core distance, as assumed in the energy reconstruction. Dashed- dotted: shower 
fluctuations and photoelectron statistics; dotted line: contribution from the error in the distance 
measurement for a constant error in core position; full line: for a core error as shown in Fig. 
4(b). 




Fig. 7: Measured mean image size as a function of core distance, for showers of fixed energy, 
using the Mkn 501 data sample (see text for details). The line shows the result of Monte Carlo 
simulations, smeared with the experimental resolution in the core distance. 
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Fig. 8: Monte-Carlo studies of the energy resolution. (a,b) Distribution in the relative difference 
between reconstructed and true energies, for two algorithms, (c) Energy resolution as a function 
of 'j-ray energy. 

measurements are correlated. Such a correlation might e.g. be caused by the variation in the 
height of the shower maximum; if this is indeed the case, it should be possible to improve 
the resolution over the values given in Fig. |H1 since the height of the shower maximum can be 
reconstructed with a resolution of about 1 radiation length (M. Ulrich et al., 1997). 



5. Effective detection area 

Different techniques are used to determine effective areas, with consistent results. In one 
approach, the effective area is determined from simulations in the usual way, imposing an 
upper limit on the distance of the core from the central telescope of 200 m. The resulting area 
(Fig.Cnta)) saturates at 7r(200 m)^ for energies above a few TeV. Alternatively, the simulations 
are used to define an energy- dependent maximum radius, such that for events within this radius 
the average trigger efficiency if large (> 80%). The resulting radius is about 100 m at 1 TeV 
and rises to 200 m at about 4 TeV; above this energy, a constant 200 m cut is applied. To 
determine spectra, in this latter case a (trigger-) efficiency correction of (10 ± 10)% is applied. 

All these estimates of effective areas rely on the ability of the simulation to reproduce 
single-telescope trigger efficiencies. These, however, can in turn be determined from system 
data. One reconstructs an n-telescope event using only a subset of n — 1 telescopes, and then 
checks if the remaining telescope has triggered. Fig. ITnT b) shows the measured trigger efficiency 
obtained for showers in two energy ranges as a function of core distance. At ~ 1 TeV, the 
trigger efficiency saturates at > 90% out to distances of 120 m. With this measured dependence, 
one can, e.g., verify that TeV showers with cores within 100 m from the central telescope will 
indeed trigger two telescopes with more than 80% average probability. 

Another way to verify the estimate the maximum radius for efficient triggering of the 
lACT system is to simply plot the distribution dn/dr"^ of events, where r is the distance to the 
central telescope. Ideally, the area density dn/dr"^ should be flat up to the (energy-dependent) 




Eij [TeV] (E.j-E^)/(E,.+E^) 



Fig. 9: (a) Scatter plot of shower energy measured with one pair of telescopes, vs. that 
measured with the other pair, for j-rays from Mkn 501 . (b ) Relative difference in the energy 
values measured by two pairs of telescopes. 
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Fig. 10: (a) (Left) Effective detection area determined from simulations, as a function of energy. 
Only events within 200 m from the central telescope are accepted (b ) (Right ) Measured single- 
telescope trigger efficiency for showers in the 1.0 - 1.2 TeV and 2-3 TeV energy range, as a 
function of the distance between the telescope and the shower core. 
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Fig. 11: Distribution dn/dr'^ of shower cores relative to the central telescope, for two energy 
ranges. 

maximum radius. Fig. ITT] shows this distribution for energies of 1 TeV to 1.5 TeV and for high 
energies above 5 TeV. The events above 5 TeV show a shght drop-off beyond radii of 150 m; 
however, the trigger efficiency averaged out to 200 m radius is still 84% (assuming that is it 
indeed 100% for small radii). At one TeV, the efficiency is flat out to 100 m. 

6. Summary 

The discussion given in this paper represent a progress report, and it is clear that we still 
have a significant way to go towards a final 'product'. On the other hand, it should be obvious 
that lACT systems provide a wealth of redundant information, both on individual air showers 
and on average shower properties, and that they promise a qualitative step forward in the 
determination of 7-ray energy spectra. The power of the HEGRA lACT system, combined 
with the large sample of about 30000 7-ray events from Mkn 501, can be used to understand 
the properties of the detection system at a level of precision which seemed impossible before. 
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